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Abstract 

We formulate the standard real-space renormalization group method in a way which takes into 
account the correlation between blocks. This is achieved in a dynamical way by means of operators 
which reflect the influence on a given block of its neighbours. We illustrate our method in the 
example of the tight-binding model in ID and 2D for various types of boundary conditions. 
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1 Introduction 



The initial success of the Real-Space Renormalization Group method applied by Wilson [Q] to the study 
of the Kondo problem rised the hope that this technique could be generalizable to more complicated 
many-body problems in Field Theory and Condensed Matter. Various groups of physicists working 
on those areas were able to combine Wilson's ideas together with Kadanoff 's concept of block thereby 
arriving to a method called Block Renormalization Group (BRG) 

From a quantitative and sometimes qualitative point of view, the BRG procedure proved to be not 
fully reliable particularly when compared with numerical approaches, such as the Quantum Montecarlo 
method, which were developed at the same time. This was one of the reasons why the BRG methods 
remained undeveloped during the 80's until 1992 were a new numerical method called the Density Matrix 
Renormalization Group (DMRG) was proposed to overcome the difficulties of the old techniques 

i- 

In order to understand the failure of the BRG method Wilson proposed in 1986 @ to study the 
problem of a single particle moving freely in a ID box. A straightforward application of the BRG 
techniques gives results which are several orders of magnitude off the exact values. According to White 
and Noack |5[ this poor performance is due to the fact that the truncation of states within each block 
in the BRG method keeps states which do not have the appropriate boundary conditions (BC). As a 
matter of fact, an isolated block has BC's which are different from those it has when it is immersed 
into the rest of the lattice. In |^] two solutions to this problem were proposed which in turn give rise 
to two different methods called Combination of Boundary Conditions (CBC) and the aforementioned 
DMRG In the CBC method one applies different BC's to each block and mix them up resulting in 
an accurate performance. However, this method is difficult to generalize for interacting systems. On 
the contrary, the DMRG method has been applied with great success to interacting systems [^, [^. In 
the DMRG method each block somehow chooses the exact BC it needs. This is achieved by putting the 
block in connection with the rest of the lattice. In practical terms, one considers a superblock where 
the block is embedded, then one finds the exact ground state of the superblock (called the target state) 
and constructs a reduced density matrix ps for the block. This density matrix pB gives the best posible 
representation of the target state on the block. Finally, one diagonalizes ps and keeps the eigenstates 
with largest eigenvalues in the truncation procedure. In this fashion, the role of the Block Hamiltonian 
Hb in the BRG method is played in the DMRG method by the density matrix, which has the virtue 
of containing the effect of the neighbourhood on a given block. 

The DMRG method has been generalized to 2D classical systems in references [jl^. In this 
framework the density matrix turns out to be intimately related to the Baxter's corner transfer matrix 
which suggests a pathway for a deeper understanding of the DMRG techniques. In reference 
fl^ we have introduced the Renormalization Quantum Group Method (qRG for short) based 



on quantum groups techniques and we mentioned that this qRG method has some analogies with the 
DMRG. 

All these works show that the DMRG is not only a powerful computational method but also a source 
of inspiration for further works concerning the RG. For these reasons, it may be worthwhile to explore 
different options or alternatives to the DMRG which may be useful in situations where the DMRG 
encounters difficulties, as in the case of 2D quantum systems. The main message of the DMRG is that 
blocks are correlated. The implementation of this idea by means of the density matrix formalism may 
be not the unique way to proceed. On the other hand, the "onion-scheme" a la Wilson adopted by the 
DMRG, while being one of the reasons of its spectacular accuracy, imposes certain limitations. 
At this stage it is not clear how fundamental are the density-matrix formalism or the onion-scheme for 
a RG method which takes into account the correlation between blocks. One can indeed combine the 
Kadanoff block method with the use of a density matrix in the process of truncation, as in reference 



14 1 . More work remains to be done to see wheather there is a real improvement of the standard 
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BRG method by combining it with the DMRG as in |14|. In this paper we want to explore another 
possibihty which is to give up both the density matrix and the onion-scheme. With this point of view 
in mind, it would seem that we should come pretty close to the standard BRG method, were it not for 
the enormous freedom hidden in a Real-Space RG method. This freedom comes from the separation 
of the Hamiltonian into an intrablock Hb and an interblock Hbb Hamiltonian. This is a source of 
ambiguities which can be sometimes mitigated with the aide of symmetry arguments, but not fully 
eliminated though. This ambiguity shows up specially for terms in the Hamiltonian acting at the 
boundaries of the block. There are no general criteria as to how to include this type of terms either into 
the intrablock or into the interblock Hamiltonians, or into both! For example, in the ID Ising model 
in a transverse field (ITF model), a choice which preserves the selfduality of the model attributes some 
self-couplings to the Hb and others to the interblock Hbb, and it yields to an exact value of the critical 
point and the critical exponent v [^], The ambiguity in the splitting of H into the sum Hb + Hbb 
thus affect deeply the truncation procedure itself, which is based on the diagonalization of Hb- Rather 
than blaming the BRG for its lack of uniqueness, we should use its freedom to allow the blocks to 
become correlated in the RG procedure. In our present approach this correlation will be taken into 
account in a "dynamical" way rather than in a "statistical" way as in the DMRG. This will be achieved 
by the introduction of interblock operators which reflect the "influence" between neighbour blocks and 
which are defined at the boundary of the block in the first step of our CBRG method. 

We have chosen to illustrate our approach the ID and 2D tight-binding models mainly for simplicity 
reasons, but we believe that our method could be applied to more complicated problems. In fact, the 
first step in this direction was already undertaken in reference where only 2 states at each stage 
of the RG-blocking were retained. This in turn allowed us to obtain the scaling law for the size 

dependence of the first-excited-state energy. 

In this paper we shall give the general mathematical structure underlying the results of reference 
|l5|| . This will allow us to retain more than two states in the RG-truncation and also to consider the 
two-dimensional tight-binding model. In this fashion, we shall recover the n? /N'^ scaling law for the 

n-th excited state of the ID model and the scaling law "^^^"^ in the 2D case. These results will then 
show that the CBRG method describes correctly the low energy behaviour of the ID and 2D Laplacian. 

This paper is organized as follows. In Sect. 2 we present the Correlated Block RG-method and apply 
it to the ID Tight-Binding model with different boundary conditions at the ends of the chain: Free-Free, 
Free-Fixed and Fixed-Fixed BC's. We correctly reproduce in each and every case the corresponding 
n^/A^^-scaling laws (iV — > oo) for the spectrum of the excited states. This is a novel result for it is 
achieved within the framework of the CBRG method in a clear and transparent fashion without having 
to resort to ad hoc mixing of different BC's states in the truncation operator involved in the CBC 
method ||5|. In Sect. 3 we extend the CBRG method to deal with two-dimensional Hamiltonians and 
apply the procedure to successfully solve the 2D tight-binding model with Free BC's. Sect. 4 is devoted 
to conclusions and prospectives. 



2 The CBRG Method: One Dimension 

The problem we want to study is the one-dimensional Tight-Binding model in an open chain with 
different boundary conditions at its ends. The Hamiltonian for this system takes the following matricial 
form. 
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where b and b' take on the values 1 (or 2) corresponding to Free (or Fixed) BC's respectively. This 
Hamiltonian is the discrete version of the Laplacian H = —9^, while the Free or Fixed BC's correspond 
in the continuum to the vanishing of the wave function (Fixed BC's) or its spatial derivative (Free 
BC's) at the ends of the chain, i.e., 



^(0) = 



11(0) 







Fixed BC 
Free BC 



(2.2) 



and similarly for b' which contains the BC at the other end of the chain. 

Hence, altogether there are 4 Hamiltonians of the type in ( |2.1| ), whose eigenstates and eigenvalues are 
the subject of our RG-techniques. 

The first step in the RG method is to divide the lattice into blocks containing Ug sites each and 
labeled with and index p {= 1, . . . , N/ug). Let us suppose for a moment that we isolate the pth-block 
from the rest of the lattice so that its dynamics, as an independent entity, is governed by a Hamiltonian 
denoted by Ap, which we may call uncorrelated block Hamiltonian. The restoration of the block back 
into the lattice involves two effects. The first one is that the BC's of the p-th block may change under 
the influence of the p+1 and p — 1 blocks. We describe this change of BC's by the action of Boundary 
Operators denoted by Bp.p±i on the pth-block. The second effect is the interaction between the pth- 
block and its neighbours p + 1 and p — 1, given by interaction Hamiltonians Cp^p±i which act on both 
p and p+1 blocks simultaneously. If the problem under consideration is translationally invariant, all 
the Hamiltonians defined above are independent of the block label p, in which case we denote them by. 



An — A 



B 

a 



p,p+i 



p,p+i 
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p,p 



Bl 



(2.3) 



The Hpree,Free Hamiltonian (2J) gives an example of this as we shall show below. Hence, for the 
time being, we shall consider the situation described by ( |2.3D and leave the more general case after 
explaining the general ideas. 

In the standard BRG method the block Hamiltonian Hb and the interblock Hamiltonian Hbb sue 
given, according to our previous definitions, by the following formulas 



Hb = A + Bl + Br 



(2.4) 



Hbb 




(2.5) 



The whole Hamiltonian is by all means the sum of Hb and Hbb for all the blocks of the chain. For a 



review on the Block RG method see [18| and chapter 11 of reference [^]. 

Next step in the RG method is to diagonalize Hb and keep its, say m {m < n^), lowest eigenstates. 
The truncation is given by a x m matrix T whose columns are precisely the components of the m 
lowest eigenstates of Hb- The renormalized Hamiltonian in the new basis is given by. 



H' = T\Hb + Hbb)T 



(2.6) 
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At first sight from Eq. (2.4) it would seem tliat we have taken into account the effect of the BC's 
on a given block. However, as the examples show, this is quite a bit illusory. On the other hand, 
the distinction among A, Bl and Bn is rather inmaterial as far as Hb is concerned, and in fact no 
distinction of this sort is made in the standard BRG formalism. Finally, let us observe that Hb and 
Hbb play rather different roles in the truncation procedure. This asymmetry has been observed as a 
source of problems by several authors in the past |17|. 

We shall mention that this asymmetry has recently been related to quantum groups in a fashion 
which has led to a new RG method called the Renormalization Quantum Group method [|l2[, [p!^]. 

Therefore, from various points of view, one is urged to make more explicit the role played by the 
BC-operators Bl and B^i in our CBRG procedure. For this purpose, we have found convenient to use 
the concept of superblock already introduced in reference [^]. We shall define a superblock as the set of 
two consecutive blocks, p and p + \ and denoted by 1). The great advantage of the superblock is 

that it allows us to materialize the distinction among A, Bl and Bf( . In fact, just as the isolation of a 
single block leads us to the definition of the Hamiltonian A, the isolation of two blocks contained in a 
superblock allows us to define B^, Bn and also C through the superblock Hamiltonian Hsb as follows. 



A + Br c 
Ct A + Bl 



HsB = \ '^;,r" .To I (2.7) 

Similarly, the Hamiltonian describing the interaction between superblocks is given by (see Fig.l) 



H. 



sB,sB 
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Br 
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Bl 



\ 



(2i 



Now instead of diagonalizing Hb in Eq. (2^), in the CBRG method we shall diagonalize Hsb in Eq. 
( |2.7| ), and afterwards keep the m = Ug lowest eigenstates in the tight-binding model. As in the standard 
BRG method, the change to the truncated basis defines the renormalized operators as follows: 



H. 



sB 



T^HsbT = A' 



(2.9) 
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sB 



T^H. 
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(2.10) 



where the matrices B'^ and C are the renormalized version of the operators A, Br, Bl and 

C, and they exhibit the same geometrical interpretation for the renormalized block as their unprimed 
partners for the original blocks. 

If we set Br = Bl = in Eqs. and (^, then after the first RG-step we get fi^ = S£ = 

and thus the previous RG-scheme coincides with the standard BRG. We may say that uncorrelated 
blocks are in a sense a fixed point of our method. However, this fixed point may be unstable, and to 
explore this possibility one has to look for non-vanishing i?-operators and their RG-evolution. 

Let us address now some examples. We shall first study the Hamiltonian ( |2.lD with Free BC's at 
the ends {b = b' = 1). Choosing = 3 for example, we see that the choice for the operators A,Br, Bl 
and C in the first step of the CBRG procedure is given by. 
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This choice is equivalent to the assumption that an isolated block satisfies Free EC's at its ends. The 
role of Bji and Bl is to join these blocks into a single chain. This is the geometrical explanation of Eqs. 
( 2.11| ). In more general cases one must have to explore which is the best choice. The generalization 



of Eqs.( |2.11 ) to blocks with more than 3 sites is obvious. In Table 1 we collect our CERG-results for 
the first 5 excited states for a chain of = 12 x 2^ = 768 sites. Comparison with exact results gives a 
good agreement. 

An important feature of our CERG method is that the r? jN'^- scaling law {N — > oo) for the 
energy of the n-excited states of a chain made up of sites, is reproduced correctly (see Fig. 2). In 
Table 2 we show the variation of the first-excited-state energy with the size N of the chain. From those 
values we can extract the corresponding 1/A^^-law which turns out to be, 

j^iflBRG)^^^ = cgLc]^' cgLc = 9-8080, (iV oo) Free-Free EC's (2.12) 

while the exact value for the proportionality constant c is c^-yr^^ox = '^^ = 9.86. This amounts to a 0.6 
% error. 

Likewise, we have enough data so as to obtain the corresponding n^/A^^-law for the whole set of 5 
excited states. Thus, the scaling law we obtain is, 

2 

^^^)(iV) = ccBRG^. ccBRG = 8.4733, {N oo) Free-Free EC's (2.13) 

which now amounts to a 7.34 % error. This is a natural fact from the worse knowledge of the highest 
excited states of the spectrum in a RG-scheme. 

We can make even more explicit the successful achievement of the 1/A^^-scaling law by leaving as 
a free adjustable parameter the exponent of \/N in addition to the proportionality constant. Let us 
denote by 9 this critical exponent. Using data from 20 to 50 steps of our CERG-method for several 
truncation of states according to our scheme 2ns i^s (namely, = 10, 13, 20) we arrive at the 
following results, 

Ef^^^^\N) = ccBRG^, (iV ^ oo) Free-Free EC's (2.14) 
For 20 — ^10 e = 1.9708 

For 26 — > 13 9 = 1.9734 (2.15) 
For 40 — > 20 9 = 1.9854 

These results clearly support the fact that we have correctly reproduced the exact value of ^ = 2 for 
the finite-size critical exponent. 

Last, but not least, as was proved in [^] our CERG method gives the exact energy of the ground state 
for every step of the RG-procedure for Free- Free EC's. 

In reference it was shown that one can reproduce easily the wave function of the excited states. 
This procedure was called reconstruction since it works "downwards" in the CERG method. The basic 
equation to be used is the reconstruction equation p^ , 

= + ij^^M (2.16) 

where ^'^''^ denotes the collection of m lowest eigenstates in the r-step of the CERG-procedure, and Ly, 
Rr are the block matrices in terms of which the truncation matrix can be written as = {Lr,Rr)- 
Our results for a chain of = 12 x 2^ = 768 sites and = 6 states kept are given in Fig. 3 where 
we have plotted the first 5 excited states and compare them with the exact wave functions. There 
are some remarkable facts regarding these figures. Firstly, the number of nodes is correctly preserved 
by our CERG wave functions. Secondly, the Free-Free type of boundary conditions are also correctly 
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reproduced at the ends of the chain. And lastly, it is worthwhile to point out that the CBRG wave 
functions "degrade gracefully" as the energy of the excited state raises in accordance with the fact that 
the lower the energy is, the more reliable are the results. 

This ends the results for the Free-Free BC's. In order to address other types of BC's we must come 
back to the case where the matrices A,Bji, Bl and C depend on each particular block. Thus, for 
example, for the Fixed- Free BC's we shall choose as the uncorrelated ^-matrix for the block located to 
the left end of the chain the following form (n^ = 3), 



A, 



I 



V 




Fixed-Free BC's (2.17) 



while the remaining matrices Ap, = 2, . . . , A^/3), will be given by Eqs.(2^), ( 2.11 ) 



For Free-Fixed BC's, it is the last ^-matrix which we have to take different from the others, namely, 

/ 1 -1 \ 

Aat/s = -1 2 -1 Free-Fixed BC's (2.18) 

As for the Fixed-Fixed BC's case, we must change the A-matrix at both ends of the chain according 
to the following prescription. 



A. 



I 



\ 





Fixed-Fixed BC's (2.19) 



Then we follow the same steps as for the Free-Free BC's, taking care that the A,Bji, Bl and C 
matrices in each CBRG-step may depend on the position of the blocks. This implies in particular that 
the embedding T-matrices may also vary from block to block. 

In Tables 3 and 4 we summarize our results for the Free-Fixed and Fixed-Fixed BC's (Fixed-Free 
BC's are equivalent to Free-Fixed BC's by parity transformation). In these tables we present our CBRG 
results for the first 6 lowest lying states for the ID tight-binding model in a chain of = 12 x 2^ = 384 
sites with mixed boundary conditions, and they are compared against the exact and standard BRG 
values. Several remarks are in order. First, we observe that the CBRG method produces a good 
agreement with the exact results and certainly much more accurate by several orders of magnitude 
than the old BRG method. Second, the CBRG method is able to reproduce the corresponding n? /N"^- 
scaling laws for the spectrum of excited states in each case of mixed BC's. Namely, 

• For Free-Fixed BC's and considering just the ground state, we have 

EiCBRG)^^^ = c^BRG^^ ^CBRG = 9-072, (A^ oo) Free-Fixed BC's (2.20) 

which amounts to ci 8 % error with, respect to the exact value of c^xact — ^ • 
As for the corresponding law for the whole spectrum, we find 

^{CBRG) ^ ^^^^^ {ri + lf ^ ccBKG = 7.6729, (AT ^ oo) Free-Fixed BC's (2.21) 
which represents a 11.5 % error with respect to the exact value of vr^. 
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• For Free-Fixed BC's and considering just the ground state, we have 

i?S™)(iV) = 4%^^^, cgJj^G = 8-35, {N oo) Fixed-Fixed BC's (2.22) 

which amounts to ci 8 % error with, respect to the exact value of c^^xact — ^ • 
As for the corresponding law for the whole spectrum, we find 

^{CBRG) ^ ^^^^^ ^"^2^^' ' ccBijG = 6.9696, (iV ^ oo) Fixed-Fixed BC's (2.23) 
which represents a 16 % error with respect to the exact value of vr^. 

We obtain bigger errors in the determination of these scahng laws as compared with the Free-Free case 
mainly because we have used less data in our fitting. Nevertherless, we find a good agreement with the 
exact results. Yet, there is another reason as to why the accuracy in the case of mixed BC's is worse, 
namely, the ground state wave function is not homogeneous in space as it is in the Free-Free case 



|15|. This makes the RG-procedure more involved and a source of extra uncertainties. 

Let us mention in passing that we are also able to make a wave function reconstruction in the mixed 

BC's cases as has been done for the Free- Free BC case. 

The outcome of all the results presented so far is that we have succeded in devising a Real-Space 
RG method capable of reproducing the correct eigenvalues and eigenstates for the tight-binding model 
as originally envisaged by Wilson, within a certain accuracy which can in principle be improved. 

Althoug the model we have employed to test our CBRG-method is a tight-binding model, there are 
some remarkable facts regarding the fixed-point structure of our CBRG-solution that we would like to 
stress. Namely, we have found that after enough number of CBRG-iterations, the matrices A, B^, Br 
and C in the Free- Free case scale nicely with the size N of the chain according to the dynamical critical 
exponent z. To be more precise, let us introduce the fixed point values of those matrices denoted by 
A* , Bl, B*^ and C* which we define as, 

A* = N-''a*, Bl = N-''bl, B*^ = N-%\, C* = N-'c*, Fixed-Point values (2.24) 

in terms of the scaled matrices o*, b^, and c*. For a block of 3 sites (n^ = 3) we find the following 
Fixed- Point structure parametrized by two constants s and t (for bigger Ug we need extra parameters). 
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s 



s t -t \ (2.25) 

s t 



with s = 1.3993 and t = 1.9581. The critical exponent z we obtain is. 



z = 0.9999 (2.26) 

which is indeed very close to the exact value z = 1 (actually, it differs in the ninth decimal digit). 
The interpretation of this Fixed-Point in the context of the CBRG method is as follows. We pointed 
out before that when the boundary operators bL^R vanish we recover the standard BRG method in 
which the blocks are not correlated. Here we find that it is the uncorrelated Hamiltonian which vanish, 
while the boundary bi^R and interaction c operators do not vanish within the scaling law. This fact 
may perhaps be interpreted by saying that in the example under study the correlation between blocks is 
more important than their selfenergy. In references |^], [10|, |11] it was shown that the DMRG method 



leads, in the thermodynamic limit, to a "product form" ansatz for the ground state wave function. In 
our case we see from Eqs.(2.24), (|2.25| ) that we also reach thermodynamical limit, which leads us to ask 
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about the nature of the ansatz for the ground state and excited states impUed by the CBRG method. 
The answer to this question will be addressed in a future publication but it suffices to say that both 
the DMRG and the CBRG methods seem to yield different ansatzs of the ground state wave function. 
In a few words, the DMRG is associated with a "vertex picture" while the CBRG is associated with a 
"string picture" . 



3 The Two-Dimensional CBRG- Algorithm 

The RG-method that we have devised in the one-dimensional problem can be generalized in a natural 
way to higher dimensions. Wc shall consider for simplicity the 2D case. First of all, we divide the square 
lattice into blocks of sites each. Each block will in turn be a square lattice with a minimum of 4 sites 
(=2x2 block). As in ID, we shall define the following Hamiltonians to carry out the CBRG-program, 

• Ap = self-energy of the p-th block isolated from the lattice. 



Bp,q = self-energy of the p-th block induced by the presence of the g-th block. 



• a 



interaction between the p-th block and the q-ih. block. 



The difference with respect to the ID case is that each block has now 4 neighbours and therefore there 
are four different B and C matrices. 

Let us consider again the Hamiltonian of a free particle moving in a 2D-box with Free BC's at the 
boundaries of the box. The 2D Hamiltonian is given again by the incidence matrix of the lattice. As 
in ID we shall choose the matrix A as the incidence matrix of the block. Thus, for example, for a 2 x 2 
block labelled as in Fig.4 we have. 



A = 



( 2 -1 
-1 2 
-1 



-1\ 
-1 
2 -1 

) 



(3.1) 



\ -1 -1 

The 4 Boundary Operators B in the same basis as in Fig.5 a have a diagramatic representation as 
shown in Fig.4 which helps us to keep track of their location in the block Hb and interblock Bbb 
Hamiltonians. Dotted line in Fig.4 means influence or correlation between neighbouring blocks. Their 
explicit matricial form is as follows. 



B\2 = S43 = Bl = 



B 



14 



B. 



23 



Bd 
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/ 1 



B2\ = S34 = Br = 



B. 



41 
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32 



Bu 



( 1 
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1/ 

0; 



(3.2) 



(3.3) 



where the labels denote the position of the neighbouring blocks and we have used the translation 
invariance of the 2D tight-binding model so that we need only to distinguish between Right and Left, 
and Up and Down. 

As for the Interaction (7-Opcrators whose diagramatic representation is depicted in Fig.4 we have the 
following matricial representation, with the same considerations as for the B-operators, 
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Clr 
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34 
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(3.4) 



Ci4 = C23 = Cdu 
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32 
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-1 \ 









(3.5) 



Thus translation invariance reduces the number of independent CBRG-matrices by a half. These 

relations are particular of the problem at hand but we must left open the posibility of having all those 

matrices different from each other in order to handle more complicated problems. 

Now that we have all the elements entering in our CBRG-method we proceed to construct the block 

HsB and interblock HsB,sB Hamiltonians out of them. To this end we have to consider a superblock 

made up of 4 blocks as shown in Fig.5 a) where the basis chosen is made explicitly. Thus, for HgB we 

have, 



HsB = 



( A + Bl + Bd Clr 

Crl A + Br + Bd 

Cud 
Cud 



(3.6) 



Cdu 


Crl 
A + Bl + Bu J 

This is a 4ns x matrix made up of Ug x Ug matrices. 

As for the interblock Hamiltonian Hsb,sB we have to distinguish between {sB, sB)-couplings of hori- 
zontal type denoted by H^'^^J^ and vertical type denoted by -ffsS sB' which read explicitly as. 





Cdu 
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Cdu 




Bu 




V Cdu 














Bu J 



(3.7) 



(3.8) 



where we have made use again of translational invariance. 

Once that we have made our choice for the decomposition of the total Hamiltonian of the 2D-tight- 
binding model into block and interblock Hamiltonians according to our CBRG-prescription, we can 
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carry on with the truncation part of the RG-method. We shall keep Ug states out of Aug states per 
superblock so that our truncation scheme may be summarized as, 

4ns (superblock) — > (new block) 

Recall that at each step of the CBRG-method we need to identify the A, Bl, Bh and C operators 
which define the truncation procedure for the next step of the method. For this purpose, firstly the 
truncation of the superblock HgB gives rise to the A' uncorrelated self-energy operator for the next 
RG-step, namely. 



HsB ( X 4ns matrix) — > A' (n^ x matrix) 



(3.9) 



To identify the rest of the operators we have to renormalize the interblock Hamiltonian which comes 
in two types, horizontal and vertical. As exemplified schematically in Fig. 6 the renormalization of the 
H^^°% Hamiltonian is given by. 



H 



ihor) 
sB,sB 



Likewise, for the H^fg^sB Hamiltonian we have (see Fig. 6), 



RL 



(3.10) 



(ver) 
sB,sB 



B'u 



(3.11) 



Now that we have identified all the operators defining the CBRG method at the new stage of the 
renormalization, we may reconstruct the new superblock Hamiltonian H'^g, which in turn has the 
same form as the original Hsb in Eq.(p.(]|) substituting all the operators by their primed versions. This 
statement can be explicitly checked by considering the set of 4 superblocks as depicted in Fig.5. Firstly, 
the new H'^^ has a contribution coming from the truncation of each of the 4 superblocks in Fig.5, each 
of them contributing with an ^'-operator as in Eq.( |3.9D . Secondly, H'^^ picks up two more contributions 
coming from the horizontal and vertical interaction between superblocks, which we denote by and 
-fft. Thus, in the CBRG-method H'^^ is renormalized as, 



( A' 



A' 



A' 



\ 



A' I 



(single superblock contribution) 



{Hi) 



( B', a. 



+ 



RL 



( B'„ 



+ 



UD 



LR 

B'r 



+ 



/ 



V 



C'du\ / 



B\ 



+ 



D 



J 



V 



B'„ a 



^R 

^LR 



RL 

B'l ) 



B' C 



DU 



UD 



B'r 



D 



and altogether we arrive at the previously stated result of Eq. 
Similarly we may proceed with the renormalized interblock Hamiltonians H'^^^Jb 



(3.12) 



O and 



(3.11) and we end up with the same form for them as the original ones. 

This ends the implementation of the CBRG-method for the 2D-tight-binding model. 
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In Table 5 we collect our CBRG results for the first 4 lowest lying states for a chain of = 
4x4x4^ = 65536 sites. Comparison with the exact results gives a good agreement. We have also 
data from truncations with blocks of = 9 and = 16 sites which enforce this statement. Moreover, 
notice that the first excited state is a doublet as in the exact solution. 

Another important result of our CBRG-method is that the (nf + n2)/N'^ scaling law for the energy 
of the (ni, n2)-excited states of a square lattice of length N is reproduced correctly. In fact, from 
data of the Ug = 16 sites truncation for the first-excited-state energy we can extract the corresponding 
1/iV^-scaling law which turns out to be, 

i^P'^^^liV) = cgJj^G]^' = 9.7365, (A^ ^ ^) D=2 Free EC's (3.13) 

while the exact value of the proportionality constant c is Cexact = vr^ = 9.86. This amounts to a 1.3 % 
error. 

Likewise, we may obtain the full {n\ + "112) /N"^ scaling law for the whole set of 15 excited states and 
we find, 

= ^cIrg ^""'^2 ^ ^CBRG = 7.9074, (iV 00) D=2 Free EC's (3.14) 

which now amounts to a 10.5 % error. 

As in the ID Free- Free case, we can determine critical scaling exponent 9 ( |2.14| ). For a truncation 
scheme 16 —s- 4 we find. 



9 = 1.99999981 D=2 (3.15) 

which clearly supports the scaling laws introduced above. Notice again (see Table 5) that our CERG 
method gives the exact (within machine precision) energy of the ground state. This is true for every 
step of the RG, as was proved in |]l5| for ID. 

We can also perform the wave function reconstruction of the excited states in the two-dimensional 
real space. This is achieved by a two-dimensional extension of the reconstruction equation ( 2.16| ). As an 
illustration of how the CERG method performs with this matter, in Fig. 7 we show a two-dimensional 
density plot of the wave function reconstruction for the siglet (l,l)-excited state in a square lattice of 
N = 65536 sites with Free EC's. In that figure we show level lines separating regions of different height 
which are colored differently (darker grey levels meaning higher height). This density plot is made by 
projecting the 3D wave function representation onto the x — y plane according to the height of the 
function. Thus, we may appreciate two maxima (black spots) and two minima (white spots) at the 
corners of the lattice, as it should be according to the exact solution. 

Recall that the exact solution for the (1,1) state has two lines of nodes, one vertical and one horizontal, 
which pass through the mid-points of each lattice side. These two node-lines are the analogue of the 
node-points of the ID excited states in Fig. 3. We have also obtained these node-lines correctly within 
our method. Thus, the qualitative real-space form of the excited-state wave functions are captured by 
the CERG procedure. 



4 Conclusions and Prospectives 

We have presented in this paper a new Real-Space Renormalization Group method called Correlated 
Elock RG (CERG) which has been very much inspired by the Density Matrix RG method introduced 
by White |5|. The basic notion developed by White in his novel treatment of Real-Spaces RG 
methods is that blocks are not isolated during the process of truncation of states and that they are 
certainly correlated. We have made a novel treatment of this concept without having to resort to a 
density matrix formalism nor to Wilson's onion-scheme for enlarging the lattice (chain). 
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The way in which we take into account correlation between blocks is a dynamical one. This means 
that we start the first step of our CBRG method defining uncorrelated block Hamiltonians A, boundary 
operators B and interaction operators C which are used to construct the block Hamiltonian Hb and 
thereby the truncation operator . This construction is reproduced at each step of the RG procedure. 
Thus, it is the system under consideration which goes choosing the correlation effects between blocks 
during the renormalization process. 

We have tested the CBRG method on the ID and 2D Tight-Binding model with different boundary 
conditions and obtained very good results for the spectrum of low lying states and its real space 
representation of their wave functions. Moreover, we have correctly reproduced the scaling laws r? jN'^ 
for the energy of the excited states in ID for Free- Free, Free-Fixed and Fixed-Fixed BC's and the scaling 
law (n^ + n2)/A^^ for Free BC's in 2D. In addition, the correct dynamical exponent z is also reproduced. 
Therefore, we believe that we have definitely solved in a transparent and systematic way the boundary 
problem for Tight-Binding models posed by Wilson in his search for a better understanding of how 
a Real-Space RG method works. 

In order to have a more extensive comprenhension of the correlations introduced in our CBRG 
method, we think it would be enlightning to make a thorough comparison between the process of 
truncation of states in our method and the one employed in the DMRG method by means of the onion- 
scheme. This issue is left for future work. Another important line of study is the extension of the 
present CBRG method to interacting systems appearing in many-body problems. In this regard, we 
are currently studying a many-body system of fermions in one dimension as another prospective use of 
our method. 
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Table captions 

Table 1 : Exact and CBRG Values of Low Lying States for the ID Tight-Binding Model for a chain 
of N = 12x2^ = 768 sites with Free-Free EC's. 

Table 2 : Exact and new CBRG values of the first excited state for the ID Tight-Binding Model 
with Free-Free BC's. 

Table 3 : Exact, Standard KG and CBRG Values of Low Lying States for the ID Tight-Binding 
Model for a chain of iV = 12 x 2^ = 384 sites with Free-Fixed BC's. 

Table 4 : Exact, Standard RG and CBRG Values of Low Lying States for the ID Tight-Binding 
Model for a chain of = 12 x 2^ = 384 sites with Fixed-Fixed BC's. 

Table 5 : Exact and CBRG Values of Low Lying States for the 2D Tight-Binding Model for a 
lattice ofAr = 4x4x46 = 65536 sites with Free BC's. 
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Figure captions 



Figure 1: Pictorical decomposition of a given Hamiltonian H into uncorrelated ^-operators, cor- 
relation Bl- Bij-operators and interaction C-operators according to the CBRG method. Bi is a su- 
perblock made up of two Li and Ri blocks. 

Figure 2 : The /N'^-\sm for the first 5 excited states of the ID Tight-Binding Model for a chain 
of A/" = 12 X 2"* sites with Free-Free EC's. This is a hiEn-him plot. 

Figure 3 : The wave function reconstruction for the first 5 excited states of the ID Tight-Binding 
Model for a chain of = 12 X 2^ = 768 sites with Free-Free EC's. We have scaled up the exact results 
by a factor of 1.23 for clarity. 

Figure 4 : Pictorical representation of uncorrelated ^-operators, correlation Bl- Br- Bjj Bd- 
operators and interaction Clr- Cj/^-operators according to the CERG method in D=2 dimensions. 

Figure 5 : a) Two-dimensional supcrblock decomposition into 4 blocks of 2 x 2 sites showing the 
choice of basis, b) Vertical and horizontal interactions between 4 neighbouring superblocks in D=2 
dimensions. 

Figure 6 : Vertical and horizontal renormalizations of neighbouring superblocks in D=2 dimen- 
sions. 

Figure 7 : Plot of the wave function reconstruction for the singlet (l,l)-excited state in the 2D 
Tight-Einding Model for a lattice of A'^ = 4 X 4 X 4^ = 65536 sites with Free EC's. 
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Energies 


Exact 


CBRG 







1.1340 X 10-^4 


El 


1.6733 X 10"^ 


1.9752 X 10"^ 


E2 


6.6932 X 10"^ 


7.6552 X 10~^ 


E3 


1.5060 X 10-"^ 


1.8041 X 10"^ 


E4 


2.6772 X 10-4 


2.9681 X 10-4 


E5 


4.1831 X 10-4 


5.1078 X 10-4 



Tabic 1: Exact and CBRG Values of Low Lying States for the ID Tight-Binding Model for a chain of 
iV = 12 X 2^ = 768 sites with Free-Free BC's. 



m 


N=12 2"^ 









12 


6.8148 X 10"^ 


6.8148 X lO"'-^ 


1 


24 


1.7110 X IC)-- 


1.7375 X Ur- 


2 


48 


4.2826 X 10"^ 


4.4694 X 10"^ 


3 


96 


1.0708 X 10"^ 


1.1515 X 10"^ 


4 


192 


2.6772 X 10-4 


2.9681 X 10-4 


5 


384 


6.6932 X 10-*^ 


7.6552 X 10-^ 


6 


768 


1.6733 X 10-^ 


1.9752 X 10-^ 




> 1 




9.8080/Ar^ 



Table 2: Exact and new CBRG values of the first excited state for the ID Tight-Binding Model with 
Free-Free BC's. 



Energies 


Exact 


Standard BRG 


CBRG 


Eq 


1.7754 X 10-^ 


1.5771 X 10-^ 


1.8409 X 10"^ 


El 


1.5043 X 10-4 


4.2679 X 10-^ 


1.6655 X 10-4 


E2 


4.1761 X 10-4 


4.2794 X 10-^ 


4.6408 X 10-4 


Es 


8.1831 X 10-4 


4.3053 X 10-^ 


9.1450 X 10-4 


E4 


1.3520 X 10"'^ 


4.3173 X 10"^ 


1.5179 X 10-^ 


E5 


2.0196 X 10-^ 


4.4288 X 10-^ 


2.2852 X 10-^ 



Table 3: Exact, Standard RG and CBRG Values of Low Lying States for the ID Tight-Binding Model 
for a chain of A/" = 12 x 2^ = 384 sites with Free-Fixed BC's. 



Eii(n'gi(^s 




Slaiidard BRG 


CBRG 




6.6585 X 10"^ 


5.8116 X 10-^ 


7.0843 X 10"^ 


El 


2.6633 X 10-4 


5.8155 X 10-^ 


2.9403 X 10-4 


E2 


5.9924 X 10-4 


5.8268 X 10-^ 


6.3690 X 10-4 


Es 


1.0653 X 10-^ 


5.8470 X 10-^ 


1.2289 X 10-^ 


E4 


1.6644 X 10-^ 


5.8717 X 10-^ 


1.7707 X 10"^ 


E5 


2.3966 X 10-^ 


5.9106 X 10-^ 


2.7311 X 10-^ 



Table 4: Exact, Standard RG and CBRG Values of Low Lying States for the ID Tight-Binding Model 
for a chain oi N = 12 x 2^ = 384 sites with Fixed-Fixed BC's. 
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Energies 


Exact 


CBRG 


Eq 





9.6114 X 10-^^ 


El 


1.5056 X 10-^ 


1.9390 X 10"'' 


E2 


1.5056 X 10-"^ 


1.9390 X 10"'' 




3.0012 X 10"^ 


3.8781 X 10"* 



Table 5: Exact and CBRG Values of Low Lying States for the 2D Tight-Binding Model for a lattice of 
N = AxAxA^ = 65536 sites with Free BC's. 
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